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I. INTRODUCTION 



The idea that the recently discovered X(3872) is a shallow molecular bound state of 
D*°D° and D*°D° mesons is extremely attractive and has motivated numerous calculations 
of X(3872) properties using effective-range theory, for a review see Ref. jl|. Going beyond 
this approximation requires including effects from dynamical pion exchange. The goal of 
this paper is to develop an effective theory of non-relativistic D mesons and pions that 
can be used to compute properties of the X(3872) systematically at low energies. Due to 
the accidental nearness of the D*-D hyperfme splitting and the pion mass, pion exchanges 
are characterized by an anomalously small scale compared to what is usually the case in 
nuclear physics [2|. We argue in this paper that, unlike in conventional nuclear physics, 
these effects can be treated using perturbation theory and compute the decay X — ► D D°ir 
to next-to-leading order (NLO) in the effective theory. 

We begin by reviewing the current experimental understanding of the X(3872). The 
X(3872) is a narrow resonance discovered by the Belle collaboration jsj in electron-positron 
collisions through the decay — > XK ± followed by the decay X —>■ J/ipn+'K^ . Its exis- 



tence has been confirmed by the CDF and D0 collaborations through its inclusive production 
in proton-antiproton collisions 0, [5] and by the Babar collaboration through the discovery 
mode B ± — > XK ± [6J|. The combined averaged mass of the X(3872) measured by these 

. n 

experiments is [7JJ 

m x = 3871.2 ± 0.5 MeV. (1) 

Note that the mass of the X(3872) is quite close to the D°D*° threshold at 3871.81 ± 0.36 
MeV {§(]. The Belle collaboration has placed an upper limit on the width of the X(3872) jjj: 

T x < 2.3 MeV (90% C.L.). (2) 

The X(3872) has also been observed in the decays X — > J/ip 7r + 7r~7r° and X — > J/ipl 
9J. The ratio of branching fractions for the three- and two-pion final states is 9j 

Br[X -> J/^tt+tt-] V ; 

Since these decays are thought to proceed through J/ipp for the J/^7r + 7r _ final state, 
and through J/ipuj for the J /if) 7r + 7r~7r final state, the ratio in Eq. ([3]) indicates a large 
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violation of isospin invariance. A near-threshold enhancement in D D ir has been observed 



in B — ► D°D°7r°K decays 10j. This is the first evidence for the decay X — > D D tt°, though 
the peak of the observed resonance in Ref. is at 3875.2 ± 0.7±?;| ± 0.8 MeV, which is 2a 
above the world-averaged X(3872) mass. The branching ratio for X — > D°D ir observed in 
Ref. [h]] is 8.8^3'g larger than the discovery mode X — > J/ijj 7r + 7r~ . The Babar collaboration 
has established Bi[X — > J/ip7i + it~} > 0.042 at 90% C.L. [111. Il2j. Various upper limits have 
been placed on the product of Br[5 ± — > XK^-] and other branching fractions of the X(3872) 
including D°D° } D + D~ Q, Xci7, Xc27, J/^ttV Q, and J/^r] Q. Upper limits have 
also been placed on the partial widths for the decay of X(3872) into e + e~ jla . \v\ and into 
77 H- 

The possible J PC quantum numbers of the X(3872) have been examined. The observation 
of X — > J/'ip'j establishes C — +. This is consistent with the shape of the tc +/ k~ invariant 
mass distributions ^, Q. Belle's angular distribution analysis of X — > J/ipTx + Tt~ favors 
J PC = 1 ++ [19]. A recent CDF analysis [2^ finds that J/ipTT + 7r~ angular distributions are 
only consistent with J PC = 1 ++ and 2 h . 

The quantum numbers J PC = 1 ++ arise if the X(3872) is a C = +, S"-wave molecular 
bound state of D°D*° + D°D*°. The possibility of a shallow molecular state is motivated 
by the proximity of the X(3872) to the D°D*° threshold and naturally explains the large 
isospin violation observed in pion decays and the dominance of the D°D°it decay mode. 
The narrow width and the non-observation of decays such as X — > xd are highly unusual 
for a conventional charmonium state above the DD threshold. From the mass in Eq. (JT]) 
and the recent measurement of the D° mass in Ref. s| one infers a binding energy 



E x = m D + m D * - m x 
= 0.6 ± 0.6 MeV. 



(4) 



This favors a bound-state interpretation of the X(3872), however, because of the large 
uncertainty, the mass alone cannot rule out a resonance or "cusp" near the D°D*° threshold 



211 ] . In this paper we will assume the A(3872) is a molecular bound state, though our 



method can be extended to the case where the A£3872) is a shallow resonance. For other 



interpretations, see Refs. 22 
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recent review can be found in Ref. 
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The interpretation as a DD* molecule is particularly predictive because the small binding 
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FIG. 1: One-pion exchange diagram for D*°D° -» D°D*° scattering. The single and double lines 
represent the spin-0 and spin-1 D mesons, respectively. The dashed line represents the 7r°. 



energy implies that the mo 



energy 39 
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ecule has universal properties that are determined by the binding 



441 ] . The small binding energy can be further exploited through 



factorization formulae for production and decay rates of the X(3872) [45 



39j and X — > D°D 0, y 4^] using the universal wavefunction 



Voloshin cal- 



culated the decays X -> D D°%° 
of the molecule. 

The main purpose of this paper is to consider the effect of 7r° exchange on the properties 
of the X(3872). Consider the one-pion exchange contribution to D*°D° — > D°D*° scattering 
depicted in Fig. [U This leads to an amplitude 



,2 ' 



(5) 



where g is the D-meson axial (transition) coupling, f n is the pion decay constant, e and e* 
are the polarization vectors of the incoming and outgoing D* mesons, respectively, and q 
is the momentum transfer. The scale fi appearing in the propagator denominator is given 
by yU 2 = A 2 — m 2 , where A is the D*-D hyperfine splitting and is the neutral pion 
mass. The hyperfine splitting, A, appears in the pion propagator because the exchanged 
pion carries energy q° ~ A as well as momentum q. Note that fi is anomalously small, 
/i Rj 45 MeV, because of the nearness of A = 142 MeV and = 135 MeV. This suggests 
that pions generate anomalously long-range effects and should be included as explicit degrees 
of freedom in the description of the molecule, if the binding energy in Eq. (J4]) is not much 
smaller than its upper limit. 

The pion interactions in the D and D* system were quantitatively analyzed using a one- 
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pion-exchange potential model by Tornqvist 47J , who actually predicted a DD* bound state 
(deuson) with a mass close to the observed A (3872). After the discovery of the A (3872), 
Swanson j^J considered a potential model that includes both a one-pion-exchange potential 
and a quark-exchange potential and found a weakly bound state in the S-wave 

jpc = 1++ 

channel. These authors worked in the isospin limit and used isospin-averaged pion masses 



and hyperfine splittings and obtained long-range Yukawa-like potentials 47|. Note that 



the effective mass term in the propagator in Eq. ([5]) has the opposite sign from what one 
typically obtains from meson exchange. This leads to a 7r°-exchange potential in position 
space which is oscillatory rather than Yukawa-like, as pointed out by Suzuki [2]. 

A central point of this paper is that the effect of it exchange can be dealt with using 
perturbation theory. Naive dimensional analysis of the relative size of two-pion and one-pion 
exchange graphs yields the ratio 



g 2 M DD ,fi 



1 

20 



1 

To 



(6) 



where Mdd* is the reduced mass of the D and D* and we have set g = 0.5 — 0.7 



This is in contrast with two-nucleon systems where a similar estimate yields 
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(7) 



where qa = 1-25 is the nucleon axial coupling and is the nucleon mass. A perturbative 
treatment of pions fails in the 3 Si channel where iteration of the spin-tensor force yields large 
corrections at next-to- next-leading order (NNLO) 53j, |54J. This is in part due to the large 
expansion parameter in Eq. (J7|) and in part due to large numerical coefficients appearing in 
the NNLO calculation. The amplitude in Eq. flS]) also gives rise to a spin-tensor force and one 
may worry that the perturbative treatment of pions will fail. However, even if large NNLO 



coefficients like those found in Ref. 



53 



54] appear in similar diagrams for the A (3872), the 



expansion parameter in Eq. ([6]) is small enough that one can reasonably expect perturbation 
theory to work. 

In this paper, we derive an effective field theory of the D°D*° and D*°D° interacting 
with neutral pions near the D°D*° threshold. This theory is very similar in structure to 



the KSW theory of NN interactions in Ref. 



51 



521 ] where a leading-order (LO) contact 



interaction is summed to all orders in perturbation theory to produce a bound state at LO 
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and pion exchange is treated perturbatively. 1 A novel feature of the effective theory for the 
X(3872) is that the hyperfine splitting of the D° and D*° is only 7 MeV above the 7r° mass 
and thus the pions are included as non-relativistic particles. In this paper we focus on the 
decay X — ► D°D°7c°. Our results are easily extended to X — > D°D 0r y. At LO our theory 



reproduces Voloshin's calculations using effective-range theory [39|, |40j. We then compute 
the NLO corrections to the decay width. These include effective-range corrections as well as 
calculable non-analytic corrections from 7r° exchange. We find that non-analytic calculable 
corrections from pion exchange are negligible and the NLO correction is dominated by 
contact interaction contributions. 

This paper is organized as follows. In Section II, we give the Lagrangian and discuss 
power counting in our theory. In Section III, we describe our calculation of the partial 
width T[X — > D°D°tt ]. In Section IV, we summarize and conclude. Appendix A describes 
how our Lagrangian is derived by integratingout the scales m n and A from heavy-hadron 



chiral perturbation theory (HH;\;PT) 57|, |58|, |59(. Appendix B gives the results of evaluating 
the individual NLO diagrams for the decay amplitude and the wavefunction renormalization. 

While this work was being completed, a related preprint [60] appeared which analyzed 
the effects of light-meson exchange on a bound state of heavy mesons near a three-meson 
threshold. This work used a scalar-meson model and calculated the entire line shape of the 
resonance to second order in the heavy-light meson coupling. Our work is complimentary 
to that of Ref. [60j in that we do not use a model but rather a Lagrangian that is directly 
relevant to the X(3872) and we go to higher order in the heavy-light meson coupling, where 
renormalization requires the introduction of higher derivative contact operators. On the 
other hand, we do not calculate the full line shape but work at the resonance peak where a 
Breit-Wigner is a suitable approximation. 



II. LAGRANGIAN AND POWER COUNTING 



The mass of the X(3872) in Eq. ([I]) is extremely close to the D°D*° threshold. Assuming 
that the X(3872) is a hadronic molecule whose constituents are a superposition of the 



1 A pionless effective theory of shallow nuclear bound states in which the leading non-derivative contact 
interaction is resummcd to all orders was first proposed in Ref. 55|. For a similar theory of the X(3872) 



see Ref. 



56] 
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D°D*° and D*°D°, the X(3872) binding energy is given by Eq. (jlj). For reasons stated 
earlier, our calculations assume positive binding energy and a molecular interpretation of 
the X(3872). The upper bound on the typical momentum of the D and D* in the bound 
state is then 7 = (2Mflo,,i?x) 1 / 2 < 48 MeV, where Mdd* is the reduced mass of the D° and 
D*°. For this binding momentum the typical velocity of the D and D* is approximately 
v D ~ (E X /2M DD ^) 1 ^ 2 < 0.02, and both the D and D* are clearly non-relativistic. We will 
use non-relativistic fields for the D and D*. 

The pion degrees of freedom are also treated non-relativistically. The maximum energy 
of the pion emitted in the decay X — > D°D V is 

2m x 

which is just 7 MeV above the 7r° mass at 134.98 MeV. The maximum pion momentum is 
approximately 44 MeV, which is comparable to both the typical D-meson momentum, po ~ 
7 ^ 48 MeV, and the momentum scale appearing in the pion-exchange graph, fi ~ 45 MeV. 
Since the velocity of the pions is v m = p^jvn^ < 0.34, a non-relativistic treatment of the pion 
fields is valid. In this respect the treatment of pions differs from ordinary chiral perturbation 



theory or the NN theory of Refs. 5l|, |52|. 



The effective Lagrangian includes the charm mesons, the anti-charm mesons, and the 
pion fields. We denote the fields that annihilate the D*°, D*°, D°, D°, and 7r° as D, D, 
D, D, and tt, respectively. To the order we are working we will not need diagrams with 
charged pions and charged D mesons so these are neglected in what follows. We construct 
an effective Lagrangian that is relevant for low-energy S'-wave DD* scattering, where the 
initial and the final states are the C = + superposition of D°D*° and D*°D°: 

\DD*) = [\D°D*°) + \D*°D )] . (9) 
v2 

An interpolating field with these quantum numbers will be used to calculate the properties 
of the V(3872). We integrate out all momentum scales much larger than the momentum 
scale set by po ~ p n ~ \i. For D mesons this corresponds to kinetic energy < 1 MeV, for 
pions the kinetic energy is < 7 MeV. The hyperfine splitting A and m n should be treated as 
large compared to the typical energy scale in the theory. We start from the Lagrangian of 



HH^PT 57J, |58|, |59j , which describes the interactions of D and D* mesons with Goldstone 



bosons, and integrate out the scales and A by rephasing fields to eliminate the large 



components of their energy. The method is similar to the rephasing used to remove the large 
mass from the energies of the fields in heavy-quark effective theory [(31]. Details are given 
in Appendix A. The effective Lagrangian is 

£ = D f yid Q + j D + ^id { 

+D f ( id + tt~~ J D + D ] lid + 77— J D + n j ( id + ^— + <5 | vr 




( — £— ] -=L= (DD f • Vtt + D f D ■ Vtt^ + h.c. 



— (DD + DDy ■ (DD + DD) 

Q 
16 



+ ^ (DD + DD) f • (D(V) 2 D + D(V) 2 D) + h.c. 

+ ^7k {BD+DD)KDM ' +Lc - + --' (10) 

where 5 = A — ~ 7 MeV. Note that /i 2 = A 2 — m 2 « 2m 7r ( i ). We use the notation 
V = V — V , and "• • • " in Eq. ( TlOl) denotes higher-order interactions. The pion decay 
constant is f n = 132 MeV with our choice of normalization. Notice that, since we are only 
interested in a C = + superposition of the D°D*° and D*°D° defined in Eq. (JHJ), contact 
interactions are written in terms of the combination of fields (DD + DD) / \/2. Because ir 
is a non-relativistic field, ir annihilates and tv creates 7r° quanta, so that the Lagrangian in 
Eq. (TIOD allows D°* -> D° + tt° and D° + tt° -> D*° transitions and forbids D*° + tt° -> D° 
and D° — > D*° + 7r°. Therefore in this effective field theory the only channels that appear are 
D* D° + D°D* and D°D 7r°. In amplitudes with external pions, we must multiply by y/2m n 
because of the normalization of the non-relativistic pion fields. In the A (3872) — > D°D ir 



decay diagrams, this will cancel the factors of 1/ v2m^ m the axial coupling and in the term 
proportional to Bi in Eq. fflUj) . 

Other channels can of course couple to the A (3872). The three-body channels D ± D°7r =F 
and D + D-7T° are above the A(3872) by only 2.8 ± 0.6 MeV and 3.0 ± 0.6 MeV, respectively. 
These channels can only appear as virtual intermediate states in A (3872) decay and self- 
energy graphs that contain at least two pion exchanges. These graphs are NNLO and 
therefore do not appear at the order we are working. 2 The D*+D" threshold lies 8.7 MeV 



This assumes that the interpolating field for the X(3872) is oc DD + DD, i.e. is constructed from neutral 
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above the X(3872), and we may integrate out these states because they lie outside the range 
of the effective theory. If kept in the theory, this intermediate state would also only appear at 
NNLO. One may worry about other nearby thresholds, especially J/ip p and J/i/jui which are 
only 1.4 ±1.1 MeV and 8.2 ±1.0 MeV, respectively, above the X(3872). The J/tp p channel 
has a much smaller energy gap than the others. However, one should take into account that 
the magnitude of the complex energy gap includes the width of the p, T p /2 = 73 MeV, 
so the J/ipp channel can be safely integrated out 4(|. A higher-precision analysis of the 
X(3872) may need to include these thresholds explicitly, especially if one wishes to describe 
the decays X — > J/ipTT + iT~ and X — > J/ip tt + 7t~7t° . We leave this to future work. 

Matching onto HH%PT yields the D°, D*°, and tt° kinetic terms as well as the axial 
D*°-D°-7T coupling. The coupling constant, g, is determined from data on the decays of 
D* mesons. The CLEO measurements of the D* + width yields g = 0.59 ± 0.07 at tree 
level 



48 



]. A NLO analysis of D* decays in Ref. 62| yields g = 0.27TJoq3- A more recent 
analysis |5(| obtains g = 0.61 at tree-level and g = 0.66 (0.53) at NLO, where the number 
outside parentheses refers to the result when virtual low-lying even-parity charmed mesons 
are included in the loop calculations and the number in parentheses refers to the result 
obtained when these states are integrated out. The uncertainty in the NLO extraction of g 
is estimated to be 20%. We will use g = 0.6 ±0.1 in this paper. 

The remaining terms in Eq. fflOl) with coefficients Co, C*2, and B\ are contact interactions 
that are not obtained from matching HH%PT but must also be included. They incorporate 
effects that come from shorter distance scales than the scale coming from 7r° exchange. 
We have only included operators needed to the order we are working. Co and C*2 mediate 
D°D*° + D°D*° scattering in the C = +, S-wave channel and have zero and two derivatives, 
respectively. B\ mediates a transition between D°D*° + D°D*° in the C — +, S- wave 
channel to a state with a D°, D°, and 7r°. 

In our power counting, po ^ Pd* ~P7r~/U~7~<3 and we calculate amplitudes in an 
expansion in powers of Q. Since the D°, D*° and 7r° are all non-relativistic, Ed ~ Ed* ~ 
E n ~ Q 2 , so the propagators of all particles are order Q~ 2 . Loop integrations are order 
Q 5 . The D*°-D°-tt axial coupling is order Q. In the exchange diagram of Fig. [H one can 



D-meson fields only. Since physical results should not depend on the choice of interpolating field, we are 
free to make this choice. 
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drop the energy dependence in the pion propagator. The factors of \ / 2m n from the vertices 
cancel the factors of 1/ (2m n ) in the momentum- dependent term in the pion propagator and 
combine with 5 to give 2m n 5 = /i 2 , reproducing the expression in Eq. (151). The pion-exchange 
amplitude is order Q° as is easily seen from Eq. ([T]). 

Only counting powers of momentum, the Feynman rules for the terms in the Lagrangian 
with coefficients Co, C2, and B\ are naively of order Q°, Q 2 and Q 1 , respectively. However, 
with this power counting the theory is perturbative and cannot produce a bound state. 
Instead we will treat Cq non-perturbatively, along the lines of Ref . 5l|, |55( , and sum diagrams 



with Co to all orders. At LO, using the power divergence subtraction (PDS) scheme one 



then finds 



c - 2n 1 mi 

ilW 7 -A PDS ' 1 ; 

where Ap^s is the dimensional-regularization parameter. 3 Taking Aprjs of order Q we find Cq 
is order Q~ l which justifies its resummation. In PDS, the coefficient C<i is order Q~ 2 as is Bi, 
as we shall see below. No other short-distance operators are needed for our NLO calculation 
of X — > D°D°tt . Feynman diagrams with C2 and B\ first contribute to X — > D°D°tt at 
NLO. 

In addition to expanding in Q, we will make one more approximation in the NLO cal- 
culation of X —>■ D°D°7T . In many cases it greatly simplifies calculations to expand in 
m^/rrio ~ 0.07. This is an approximation we will perform when evaluating loop diagrams. 
It is not systematized in our power counting scheme. 

As emphasized earlier, the perturbative character of pion exchange depends on the small- 
ness of the parameter appearing in Eq. ([6]). Our effective theory can be used even if dimen- 
sionless parameters conspire to render pion exchange non-perturbative, but in this case 
one-pion exchange would have to be resummed as done in the NN system 63]. 



III. DECAY RATE FOR X(3872) -> D D tt° 

Here we describe our method for calculating the width of the A(3872) resonance. We 
consider the following two-point function of interpolating fields X 1 = (D° D * 1 + D° D * 1 ) / y/2 

3 The dimensional-regularization parameter is usually denoted /i but we use a different symbol here to avoid 
confusion with the scale appearing in pion exchange. 
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for the X(3872) with the spin index i: 

G{E)5 ij = j d 4 xe-^(0\T[X\x)X^0)}\0) = E + ^/*l T/2 + - > ( 12 ) 

where Ex is the binding energy of the X(3872) and the ellipsis represents terms that are 
less important in the resonance region, E + Ex ~ T. We can define a function £(£7), where 
— z£(— represents the Co-irreducible graphs contributing to G(E). Our definition of 
E(-B) is similar to the function £ defined in Appendix A of Ref. 64j. In terms of S(-E'), 
is 

-iE(.E) 



G(E) 



C £(£) 

-iRe£(£) + Im£(£) 



(13) 



1 + C Re £(£) + z C Im £(£) ' 

Since the real part of the denominator must vanish at E — —Ex, we have 1 + 
Cq Re £(— Ex) = 0, and expanding about E = —Ex we obtain for G(E) 

i{l/C — (E + E x )Re T!{-E x )) + ImE(-£ x ) 



G(E) 



C (E + E x )Re Z'{-E x ) + i C ImZ(-E x ) 

i i 



C 2 (E + E x )Re Z'(-Ex) + % C 2 Im E(-E x ) C 
where £' = dY,/dE. From Eq. fTB"|) . we immediately see that 



(14) 



1 2Im£(-£ x ) 
( } ~ C 2 Re£'(-£ x ) ' Re £'(-£*) " (15) 

The function 2 Im £ corresponds to the square of the decay diagrams. It is interesting to 
compare the result of evaluating the loop diagrams and taking the real part with direct 
evaluation of the decay diagrams. 

Consider for example the evaluation of the two-loop diagram in Fig. [6^ in Appendix B. 
The result of evaluating the graph is 

2 i /A \ 8-2D 

g 2 1 /ApDS X 



F «- & > = - (16) 



X 



d D q f d D l 1 1 



(2tt) b J (2u) D q + E x /2 - q 2 /(2m D *) + te -q + E x /2 - q 2 /(2m D ) + it 
1 1 



x 



k + E x /2 - l 2 /(2m D *) + te -l + E x /2 - l 2 /(2m D ) + te 

(g + Mq + 1) j 

qo + k-(q + l) 2 /(2m w ) +S + ie' 

11 



We perform the energy integrals by contour integration, taking the poles of the D-meson 
propagators. This yields 

* 8-2D 

Fig.Ek) 

*JZ \ Z / 

d D ~H 



.9 2 


( A r) 








dP~ x q 


(2tt) d - 1 



(q + l)i(q + l) 3 

X 



2m n (E x - q 2 /(2m D ) - l 2 /(2m D )) - (q + I) 2 + /i 2 + ie 



,2 / A \ 8-2D 



■ 9 /ofl/T \2 / ^PDS 



2/| V 2 

/" rf 15 - 1 / 1 1 



(271) - 1 J (2tt) d ~ 1 q 2 + -f 2 - ie I 2 + 7 2 - ie 

x (g + 0t(g + Qj (17) 

-m^ 2 /M DD , +q 2 /m D + l 2 /m D ))-(q + l) 2 + /j, 2 + ie' 1 J 

The first two propagators that come from the D* mesons clearly scale as Q~ 2 . The last 
two terms in the pion propagator denominator, — (q + I) 2 + /i 2 , scale as Q 2 while the other 
terms scale as (m n /m£,)Q 2 . Since m^/m^, ~ 0.07 is comparable to our expansion parameter 
in Eq. ([6]), these terms can be systematically dropped. The neglected terms come from 
the pion kinetic energy, and in dropping them we are treating the pions in the potential 
approximation [65] . The final answer is then 

n 2 /A \ 8-2D 

Fig. Eh) = -^(2ilW 2 ( ApDS 



x 



2/r uu 1 \ 2 

d°- l q f d°-H 1 1 (q + l)i(q + l)j 



(2tt) d ~ 1 J (2tt) d - 1 q 2 + -f 2 - ie I 2 + 7 2 - ie (q + I) 2 - /i 2 - ie 



g 2 Sij ( M DD N 1 



2/2 3 V 2vr 



2 2/1,1,1 / ^-PDS 



(Apds-7) +H — + - + log 



4e 2 °V2 7 -i/i 



.(18) 



where 1 / (4e) = 1/ (4e) + (In 7r — 7e)/2. We have used three-dimensional rotational invariance 
to replace (q + l)i(q + l)j with (q + l) 2 5ij/3 and use the PDS scheme to evaluate the remaining 
scalar integrals. 

There is one instance when dropping m^/m^ corrections is not appropriate. To see this 
consider evaluating the real part of Fig. by evaluating the cut diagram. The cut runs 
through the D-meson and pion propagators. In the cut diagrams, these propagators are 
replaced with 5-functions. For the D-meson propagators integrating over the (^-functions is 
equivalent to taking the pole using contour integration. So the cut diagram is obtained from 
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Eq. (jTTI) simply by replacing the pion propagator with the corresponding ^-function. Doing 
this and making the substitutions q — > p D and I — > pd so that q + I = Pd + Pd = ~Pn, one 
obtains 

9 2 (0A/r ^2 f d 3 p D d 3 p D 

■\IM DD * 



2/r ' j w 

x |e- PA* 2 T 2 5 1 2 ^ (V " J* ~ - — P 2 , - —4) • (19) 

This clearly reproduces the interference term in Voloshin's effective-range calculation of 
X — > D°D 7r° 39(. The 5-function in Eq. ( 1T91) imposes the constraint on the phase space due 



to energy conservation. Dropping m^/m^-suppressed terms in the 5-function corresponds to 
neglecting the final-state D-meson's kinetic energy and would leave the integrals over their 
momentum unconstrained. Clearly this is not a good approximation. Physically, it is also 
clear that the on-shell propagating pion in the final state cannot be treated in the potential 
approximation. 

Therefore, in evaluating ImS(- Ex) we will calculate the decay amplitudes for the dia- 
grams and integrate over the physical three-body phase. In diagrams with virtual pions, we 
drop the kinetic energy so the pions are potential. 4 In the virtual diagrams this approxi- 
mation is valid up to 0{m lx /m£ ) ) corrections. Since our expansion parameter is expected to 
be 0.05 — 0.1, making this approximation in the virtual NLO graphs induces an error of the 
same size as the NNLO correction. 

The LO decay diagram is shown in Fig. [2j The D* propagator scales as l/Q 2 and the 
axial coupling scales as Q so the LO diagram is order Q~ x . We show only one diagram, but 
there are two channels related by C-conjugation that are implied. It is straightforward to 
evaluate these diagrams and obtain 

. g M DD * _ 

»T 2 T 2 Pt ■ ex + (Pd -> Pd) , (20) 
J-K Pd + 7 

The LO contribution to the wavefunction diagram is shown in Fig. [3j The graph is 0{Q) 
and therefore ReS'(— Ex) is 0(Q^ 1 ). The result of evaluating this graph and taking the 
derivative is 

K^oA- (2D 



For further discussion on the role of recoil corrections, see Ref. [6 
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FIG. 2: LO diagram for decay rate. 




FIG. 3: LO diagram for calculating wavefunction renormalization. 



The LO decay diagram in Fig. [2] is 0(Q~ l ) so the leading contribution to ImS(-Ex) from 
Fig.[2]is 0(Q~ 2 ). Dividing by the LO wavefunction renormalization which is 0{Q^ V ) one sees 
that the leading contribution to the decay rate is OiQ^ 1 ). The result reproduces Voloshin's 



calculation of X 



-> D°D°ir 
dT LO 



39|]: 



dp 2 D dp 2 D 



2tt7(p; • e x ) 



1 



1 



.Pd + 1 Pn + T. 



(22) 



The NLO corrections to the decay rate are suppressed by one power of Q. They come 
from graphs, shown in Figs. H] and [5j with one additional pion exchange or one insertion 
of C2 and B\. These coefficients scale as Q~ 2 . NLO contributions to the wavefunction 
renormalization are down by one power of Q as well. These contributions are given by the 
two-loop self-energy diagrams involving pion exchange or an insertion of C2 shown in Fig. [HI 

The results for individual diagrams are given in Appendix B. The final expression for the 
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FIG. 4: NLO diagrams for the decay rate involving pion exchange. 



a) 




b) 




FIG. 5: NLO diagrams for the decay rate which involve contact interaction. 



NLO differential rate is 
rfr N T,n dlYn ( ^ g 2 M 1 



dp 2 D dp 2 D 



dp 2 D dp 2 j 



1 + 



DD* 



7 ^4 7 2 -/i 2 \ ,M DD « 7 ( 7 

+ tv2l A PDSj 



gry ( gM DD > 



A 



PDSj 



# 4 M D£ ,* 7 



(Ptt ■ ex 



6tt/ 2 \4 7 2 + fi 2 

C 2 (Apds) + ^i(Apds) J (Apds - 7) (Pn ■ e x ) 2 

1 



7T 



(23) 



.Pd + 7 2 



Pn + 7 2 



64^3/4 

5f 4 M DD . 7 \Reh 2 (p D 



Rehi(p D ) Reh^pD 

P 2 D + 1 2 



P 2 n + 7 2 



Pd + I 2 



+ 



P 2 r> + 7 2 



647T3/4 



Pn-exPD- e X Pn-PD+ (j>D ~> P£>) 



+ 



+ 7 + 7 



Pd + 7 2 

The functions /ii(p) and ^(p) are given in Appendix B. The first line in Eq. 
multiplicative correction to the LO decay rate. Note that in the absence of pions 

2tt r 1 



C 2 (A 



PDSJ 



is a 



(24) 



M DD . 2 (Ap DS - 7 ) ; 

where r is the effective range. The term proportional to C 2 in the first line of Eq. fl23|) 
reproduces the expected correction from the effective-range theory, in which the leading 
correction involving r comes from the modification of the normalization of the wavefunction: 

^ ER (r) = 



7 



-jr 



— 7 r ) r 



(25) 
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CP 

c) 



FIG. 6: NLO diagrams for calculating wavefunction renormalization. 

The second line in Eq. ( 123]) is the interference between a short-distance local coupling of the 
X to the D°D°tt state and the LO amplitude. Note the coefficient of this term scale as 
1/(Apds — 7) and disappears if one takes Aprjs — * 00, confirming the short-distance nature of 
the contribution. The final terms are non-analytic corrections due to pion exchange. These 
contributions turn out to give a very small (~ 1%) contribution to the decay rate, so the 
NLO correction is entirely dominated by the contact interaction contributions. 

We will parametrize C2 according to Eq. (|24|) . where tq is to be interpreted as the short- 
distance contribution to the effective range. Since we have integrated out the scales and 
A, it is reasonable to take r$ ~ (100 MeV) -1 . We will parametrize 

<j}[lW *C 2 (A Pm ) + B 1 (A PDS )) (A PDS - 7) = 777^7 TTT, . ( 26 ) 



U zv ™' ' " (100 MeV) 3 

where rj is a dimensionless parameter we expect to be of order unity. Fig. [7] shows the partial 
width F[X — > D°D°'K ] as a function of the binding energy. The central solid line is the LO 
result. We use the central value for the tree- level extraction of the D-meson axial coupling, 
g = 0.6. The band in Fig. [7J shows the NLO rate with the parameters r and rj varied 
between 

< r < 1 , -Kn<l. (27) 

- - 100 MeV ' - ' - K J 

As stated earlier the non-analytic calculable corrections from pion exchange in Eq. (123]) give 
negligible corrections. The band is dominated entirely by the contact interaction contribu- 
tions. Measurements of the X mass and partial decay width into D°D°tt° can naturally be 
explained within a molecular picture if the corresponding point in Fig. [7J falls within, or - 
due to higher orders — close to, this band. Values far outside the band can be acommodated 
only if short-range parameters or higher-order effects are anomalously large. In either case 
the appeal of our framework would be strongly diminished. 
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FIG. 7: Decay rate for X — > D D°ir as a function of Ex- We use g = 0.6. The central solid 
line corresponds to the LO prediction. The band is the result of the NLO calculation when the 
parameters r$ and r\ are varied in the ranges < r < (100 MeV)" 1 and -1 < rj < 1. 



IV. SUMMARY 



In this paper we have developed an effective field theory of non-relativistic pions and D 
mesons that can be used to describe the properties of the X(3872), assuming it is a weakly 
bound state of D°D*° and D*°D° with anomalously small binding energy. Because of an 
accidental cancellation between the D-meson hyperfine splitting and the mass of the tt°, pion 
exchange is characterized by a smaller scale than is typically the case in nuclear physics. 
This relatively small scale and the small axial coupling in the D-meson system (compared 
to the nucleon's axial coupling) combine to make the corrections from 7r°-meson exchange 
amenable to perturbation theory. This justifies the application of a theory similar to that 
proposed by Kaplan, Savage, and Wise for low-energy NN interactions 5jJ, l52| . in which a 
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leading-order contact interaction is resummed to all orders to produce the bound state, and 
pion exchange and higher- derivative contact interactions are treated within perturbation 
theory. 

This theory reproduces at leading order the calculation of T[X — ► D°D°tt ] by Voloshin 



39( | which exploits the universal behavior of the DD* wavef unction in limit of small binding 



energy. Effective-range corrections as well as other corrections from short- distance scales are 
encoded in higher-dimension contact operators in the theory. These corrections turn out to 
completely dominate non-analytic calculable corrections from 7r° exchange. Varying these 
coefficients within ranges determined by naturalness allows us to estimate the size of correc- 
tions to the leading-order calculations of Voloshin. While it is somewhat disappointing that 
the non-analytic calculable corrections from 7r° exchange are so small that an experimental 
test of this aspect of the theory seems unlikely in the foreseeable future, the smallness of 
these corrections confirms one of the main points of this work, namely that pion exchange 
can be dealt with using perturbation theory. 

A naive estimate of the size of the NNLO corrections based on the expansion parameter 
in Eq. ([6]) is 1% or smaller. It is important to remember that in conventional nuclear 
physics, large corrections come from graphs with two or more pion exchanges in the 3 Si 
channel, which first arise at NNLO. The two-pion exchange graphs at NNLO graphs come 
with large coefficients, ~ 5, which ruin the perturbative expansion of KSW for two-nucleon 
systems [5^ |. In our case similar size coefficients in two-pion exchange graphs should not 
ruin perturbation theory since even with a large coefficient ~ 5, they would only be expected 
to be 5% or smaller. It would be interesting to perform the NNLO calculation to check this. 
A NNLO correction of 5% would dominate the non-analytic NLO contribution but would be 
smaller than the uncertainty in the contact interaction contribution, indicating convergence 
of the expansion. In the unlikely case that pion exchange is non-perturbative, it can be 
resummed as done in nuclear physics [631 ] . 

It is straightforward to extend the analysis of this paper to other V(3872) decay and 
production processes, such as X — > D°D ^ or X — > J/ipp* — > J/ipir+n^. Coupling to J/ip p 
and J/ipu channels can be incorporated by including these degrees of freedom explicitly 
in the theory and coupling them to D°D*° + D*°D° via contact interactions. It would be 
interesting to calculate tt° exchange to other decays of the V(3872) to see if these correc- 
tions lead to any interesting observable effects. It would also be interesting to use data or 
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theoretical calculations to fix some of the counterterms appearing in the theory so as render 
calculations in this paper more predictive. 
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APPENDIX A: DERIVING THE EFT LAGRANGIAN FROM HHxPT 



Heavy-hadron chiral perturbation theory (HH%PT ) 57|, |58|, |59( can be used to derive 
the low-energy effective Lagrangian for D mesons, D* mesons, and pions relevant to the 
AT (3872). We begin with the two-component HH;\RT Lagrangian introduced in Ref. 67j |. 



£ = Tr [H\iD )H] - g Tr [H ] Hct ■ A] + — Tr [rftrHa] 



(Al) 



where er are the Pauli matrices, A = — m#, H = D-cr + D, and A = — Vtt/ f n + 0(7T 3 ). 
Here D is a heavy vector field, D is a heavy pseudoscalar field, and tt is the pion field, 



7T 



71 



7T 



(A2) 



Evaluating the traces in Eq. ( 1A1I) we obtain 

C = 2D f (iD - j^j D + 2D^ (iD + ^ D 

-2g (D f • AD + D j B ■ A) - 2ig D f • D x A . 



(A3) 
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Since we wish to describe a bound state of two heavy mesons the power counting of HH^PT 
in powers of 1/ tuh is inappropriate. Instead we need to power count in the relative velocity 
v <C 1 of the heavy mesons. The kinetic energy which is sub-leading in is leading in 

v , and as a consequence we must include the kinetic term in our Lagrangian, 

J V 2 a\ J V 2 3A\ 

—2g (D" 1 " • AD + D^D ■A)-2ijD t -DxA. (A4) 

We now rescale the heavy-meson fields 

{ D iDH i, e ^{ AD}) (A5) 

which gives 

£ = Dt (^ + 2^r- A ) D + ct (> D ° + £) D 

—g (D" 1 " • AD + D^D • A) — D" 1 " ■ D x A . (A6) 

Since we are only interested in those terms involving D*°, D , and tt° we keep only these 
fields in the Lagrangian, 

£ = Dt ( !B ° + 2£:- A ) D + B, (; D <> + &}) D 

+— £— (d* • Vtt /) + D^D ■ Vtt°) - i— £— D f £> x Vtt . (A7) 
v2/,r ^ ' v2/tt 

Next the kinetic term for the pion is derived from the chiral Lagrangian 

= ^ Tr [d^d^] + lfB Ti [M (S + St)] (A8) 

4 

1 1 

= T^d^d^n + -m 2 (7r ) 2 + self-interactions (A9) 
= I^o (_^2 _ m y + . . . ) (A10) 

where £ = exp(2i7r// 7r ), A^ = diag(m u ,md) is the quark-mass matrix and Bq is a constant. 
The pion self-interaction terms are not needed at the order we are working so they are 
dropped, and we add the pion kinetic term in the last line above to Eq. (|A7[) to obtain our 
Lagrangian. 

However, this Lagrangian still includes the large scales m n and A, which must be inte- 
grated out of the theory. Since we are interested in a non-relativistic theory of pions we are 

20 



justified in splitting the pion fields into creation and annihilation operators 7r° = tt + it'. In 
addition we rescale the meson fields to make the large scales explicit: 

tt = -=L= e~ m ^ Ti , 7r j = —^e im ^7i j , D^e" im ^D. (All) 

The pion kinetic term can then be expanded 

At = ^° {-91 + V 2 - ml) vr 

= — ^ (2im w d - d 2 Q + V 2 ) ti + tt (-2mv<9 - <9 2 + V 2 ) tt* 

+e- 2im ^7r ^im^Sb - d 2 Q + V 2 ) tt + e^V (-2im*Sb - dl + V 2 ) tt*} 
/ V 2 \ 

= 7r^ i(9o H ) 7r + higher-order relativistic corrections . (A12) 

The terms in the third line include a large phase factor and as a consequence can be inte- 
grated out. In addition to modifying the pion propagator the field redefinition in Eq. flAl lj) 
modifies the kinetic term for D, 

D f (ido + - A ] D - Dt (id + - 5) D, (A13) 

where 5 = A — ~ 7 MeV. Note that after the field rescaling the last term in Eq. (IA7|) 
contains a phase factor e~* m,rt , and can be dropped. 

Finally, we obtain the the first three lines of the Lagrangian given in Eq. (fTUI) by another 
rephasing of the D and n fields, 

D-»e~ wt D, 7r — ► e~^*7r. (A14) 

This just shifts the residual mass from the D kinetic term to the n kinetic term. This 
last step is not essential, however it is convenient. The remaining terms in Eq. (fTU|) are 
short-distance interactions allowed by power counting and the symmetries of the theory. 



APPENDIX B: NLO DIAGRAMS 



The NLO decay diagrams involving pion exchange are shown in Fig. HI The result for 
Fig. Hk) is 



a 3 M 2 

y ivi dd* 



-Ap D s - hi(p D ) ) p 1T -e x + h 2 (p D ) Pd • e x Pd ■ P-k 



+ (pi>- Pd),(B1) 
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where the functions h\(p) and /12(f) are given by: 

hx(p) = / dx\J—p 2 x 2 + (p 2 + 7 2 + /i 2 )x — /i 2 — ie 
Jo 

p 2 - 7 2 - u 2 (p 2 + 7 2 + u 2 ) 2 - 4p 2 u 2 , 207 
-7 H — tan 



^(p) = dx 



4p 2 8p 3 —p 2 + 7 2 + /i 2 

• 2^ 2n . (p 2 + 7 2 + /i 2 ) 2 -4p 2 /i 2 7 2 + (^~P) 2 , R9 ^ 

~V (P +7 IP ln 7 2 + (^ + P ) 2 (B2) 

1 x 2 



v/— • p 2 x 2 + (p 2 + 7 2 + /i 2 )a; — fi 2 — ie 
5p 2 + 37 2 + 3/i 2 3(p 2 + 7 2 + ;u 2 ) 2 - 4p 2 /i 2 _ x 2p7 



7 H ;tt tan 



4p 4 ' ' 8p 5 -p 2 + 7 2 + /i 2 

,3/i 2 2 2 , 3(p 2 + 7 2 + /i 2 ) 2 -4pV 1 7 2 + - P) 2 
V (P +7 + 1^ N 2 + (/* + p) 2 ' 



and the result for Fig. [4b ) is 

g 3 M 2 DD , fi 3 



10 ta /■ 2 2^2 • ex + (pd -*Pfl)- ( B4 ) 
12tt/3 (^ + 7 2)2 

Note that Fig. Hb) includes, as a subgraph, the one-loop D* self-energy contribution. In the 
PDS scheme the self-energy graph has a linear divergence which gives an additive renormal- 
ization to the residual mass term of the D*. At tree level, we performed a field redefinition 
which moved the residual mass term from the kinetic term of the D* to the kinetic term 
of the pion through a field redefinition. In order that loop corrections do not reintroduce a 
residual mass for the D*, we introduce a counterterm —S^D^D + D), which is defined 
to cancel the residual mass term at each order in perturbation theory. At one-loop order, 
<5ct = gV 2 Apds /247r/ 2 . This linearly divergent contribution to the self-energy also appears 
in Fig. [6b) and is canceled by an insertion of the residual mass counterterm in a DD* bubble 
(not shown in the figure). 

The NLO diagrams with the counterterms C2 and B\ are shown in Fig. EJ The result for 
Fig. Eh) is 



( gM 2 D ,(A PDS -7) p^-7 2 
4vr/^ p 2 D + 7 2 



-2C 2 (Apds) — £ 2 I " e * + kPd -> Pd) » ( B5 ) 



and the result for Fig. Eb) is 



■ rj / a \ MlXD* (ApDS ~ 7) - 

- iB x (Apds) = Ptt ■ ex • (B6) 

Z7T 



22 



Finally, we show the NLO wave function renormalization diagrams in Fig. [61 The NLO 
contribution to ReS'(— Ex) from the graphs in Fig. [6] is 



Apps - 7 

7 



+ 



2fi 2 
47 2 + {i 2 



-C 2 (A 



PDSj 



M|, D »(7-Ap DS )(27-Ap D s) 
2tt 2 



(B7) 
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